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Abstract. For a homogeneous random walk in the quarter plane with nearest-neighbor 
transitions, starting from some state (io,Jo), we study the event that the walk reaches 
the vertical axis, before reaching the horizontal axis. We derive a certain integral 
representation for the probability of this event, and an asymptotic expression for the 
case when io becomes large, a situation in which the event becomes highly unlikely. 
The integral representation follows from the solution of a boundary value problem and 
is involves a conformal gluing function. The asymptotic expression follows from the 
asymptotic evaluation of this integral. Our results find applications in a model for 
nucleosome shifting, the voter model and the asymmetric exclusion process. 



1. Introduction 

Consider homogeneous random walks in the quarter plane with nearest-neighbor 
transitions. For such random walks, starting from some state (io,jo), we study the event 
of reaching the vertical axis, before reaching the horizontal axis. We derive an integral 
representation for the probability of this event, and an asymptotic expression for the case 
when iq becomes large, a situation in which the event becomes highly unlikely. We use 
the classical method of solving for the generating function via functional equations and 
boundary value problems. 

Our primary motivation is the work of Opheusden and Redig [16] on the following one- 
dimensional particle system. Consider three particles in Z and let r]£ (n) denote the position 
of particle I after n steps, with initial positions 771(0) < 772(0) < 773(0). The particles each 
get a weight and are then equipped with the following dynamics. At each time step, one 
of the particles is selected with probabilities proportional to their weights. The chosen 
particle is then moved to either the left or the right, with equal probability. Denote by 
X(n) = 772(71) — 771(77) and Y{n) = 773(77) — 772(77) the pairwise distances. The discrete-time 
Markov chain (X(n), V(n)) n6 z + , with Z + = {0, 1, . . .}, then clearly is a random walk in 
the quarter plane In particular, with particle two having weight v and particles one 
and three weights A, one obtains the walk in Figure 1(a). For this walk Opheusden and 
Redig [16] studied the event of the Markov chain, starting from (X(0),Y(0)) = (14,1), 
reaching the vertical axis, before the horizontal axis. This event plays an important role in 
studying a nucleosome shifting with respect to the DNA sequence. In [16] an asymptotic 
expression was derived; see (2.3). 

Another application of this work is the one-dimensional voter model. This lattice-based 
interacting particle system is used to model the spread of an opinion through a static 
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FIGURE 1. Walks considered in [16], [2] and [9], respectively 



population via nearest-neighbor interactions, and finds application in modeling competing 
species. It is a discrete-time process on {0, 1} Z , where each site in Z is labeled either 
or 1. Two adjacent sites (a pair) are called an unlike pair when the labels are 01 or 10. 
The voter model then has the following dynamics. At each time step the model selects 
uniformly at random from amongst all unlike pairs. The chosen pair is flipped to either 00 
or 11, with equal chance of each. Since our version of the voter model lives on the infinite 
lattice, the ground state is the so-called Heaviside configuration 

. . . 111000. . . , 

and the initial configuration is assumed to be one with a finite number of unlike pairs. For 
example, with A the number of finite blocks of zeros (or ones), a configuration looks like 

K\ ai K2 <T2 KJV <*N 

. . . in o^rm^Trnii . . . ooo^rmooo 

with Kg (resp. ag) the size of the £th block of zeros (resp. ones). It is clear that A is a 
non-increasing function of time, since blocks will merge as time progresses. In fact, the 
Heaviside configuration (A = 0) is an absorbing state. Therefore, for the voter model, 
a crucial characteristic is the hitting time r of the Heaviside configuration. In [2] it is 
shown that E[t 3 / 2_ie ] < oo and E[r 3 / 2+e ] = oo, for any e > and any initial configuration. 
In proving the latter fact, it suffices to consider the case N = 1, because A is non- 
increasing and hence the process always has to pass before absorption through A = 1. 
Therefore, in [2] the process («i(ra), 0"i(n)) n gz + = (X(n), Y (n)) n6 ^ + is considered, which 
is clearly a random walk in the quarter plane that is absorbed when it reaches the boundary 
{(0,0)} U {(i,0) : i 5s 1} U {(0, j) : j ^ 1}. Define Z* + = {1,2,...}. With 

p id = P [(X(n + 1), Y(n + 1)) = (X(n),Y(n)) + (i,j)\(X(n),Y(n)) G Zf] , 

the dynamics of the voter model is described by pi.o = Pi-i = Po,-i = P-i,0 = P—i,i = 
Po.i = 1/6, see Figure 1(b). This random walk thus plays an important role in the voter 
model. It describes the situation in which the final two remaining groups try to impose 
their opinions on each other. This situation is in many cases rather persistent, particularly 
when both groups are of considerable size; see [3] for some simulation results that support 
this fact. We are interested in the situation in which one of the two groups forms a clear 
minority, and nevertheless, wins the battle with the other much larger group. It is clear that 
this is a large deviations event, and it is for this event that we obtain precise asymptotics. 
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A third application of the results in this paper is the phenomenon of spontaneous 
symmetry breaking in the asymmetric exclusion process. Godreche et al. [9] show that 
in some limiting regime, this phenomenon can be formulated as the hitting probability 
of the random walk with the transitions as in Figure 1(c). Contrary to the first two 
applications, this random walk clearly has a negative drift. 

Hence, in all three applications, we are interested in random walks reaching the vertical 
axis, before the horizontal axis, in situations where reaching first the horizontal axis is much 
more likely. In the next section we present our results for the probability of this event, for 
both zero-drift and negative-drift random walks. We shall also discuss the consequences 
for the three applications. 

2. Main results 

Denote by (X,Y) = (X(n),Y{n)) n ^z + a random walk in the quarter plane and let 
P(j J0 ) [d?] be the probability of event S conditional on (X(0),Y(0)) = (io,jo). Throughout 
we shall make the following assumption: 

(HI) The walk is homogeneous inside of the quarter plane, with transition probabilities 
{pi,j}—l^ij^l to the eight nearest neighbors. 
Denote the horizontal and vertical axes by 

n = {(i,0) :i>0}, V = {(0,j):j^0}, 

and define H* = T~L \ {(0, 0)} and V* = V \ {(0, 0)}. The principal object of study in this 
paper is the probability 

(2.1) ¥ {iQ>jo) [(X,Y) hits V before %*}, 

for which we derive an exact expression, as well as an asymptotic expression for the large- 
deviations case iq — > oo. 

In this paper we shall for the most part restrict to random walks (X, Y) that, besides 
(HI), satisfy the following assumptions: 

(H2) In the list pi > i,pi > o,Pi ) -i,.Po,-i)P-i,-i>.P-i,0).P-i,i)R>,i) there are no three consec- 
utive zeros; 
(H3) pi,i +p-i-i +pi -i < 1; 

(H4) The drifts are non-positive: Y2—i<i jKi^Pid ^ anc ^ S-i<j j<iJPi.j ^ 0- 
Assumption (H4) guarantees that the random walk will hit one of the boundaries with 
probability one. With assumption (HI) and (H3) we can use the general framework for 
random walks in the quarter plane developed by Fayolle et al. [6] (see Subsection 3.4 for 
more details). Assumption (H2) excludes degenerate random walks, which can typically 
be analyzed using easier methods. 
Here is our first main result. 

Theorem 1. Let (X, Y) be a random walk satisfying (H1)-(H4). // 

( 2 - 2 ) E-Kij^i m,j = and E-Ki^iiftj =0, 

there exists a constant A S (0, oo) such that 

P {iodo) [(X,Y) hits V before %*} ~ A—, i oo. 
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In Section 4 we present an explicit expression for the constant A in Theorem 1. For the 
model in Figure 1(a) of nucleosome shifting, an inspection of Theorem 1 reveals 



1 — — ■ 

(2.3) P (ioJo) [(X,Y) hits V before H*] ~ -* P^r~, k -> oo. 

k ' arccos(^) z 

Opheusden and Redig [16] were able to derive (2.3) using the following approach. First, 
define a generating function of which the probabilities in (2.1) are the coefficients. Then 
show that this generating function satisfies a certain functional equation. This is the 
functional equation that is archetypal of random walks in the quarter plane, see (3.11). 
The functional equation defines a characteristic curve, and by considering its tangent points 
one can determine the dominant singularity of the generating function. The nature of this 
dominant singularity then gives the asymptotic decay, in this case 0{1/iq). Hence, an 
asymptotic estimate for (2.1) is derived by studying a functional equation without having 
to solve it. The only drawback of this approach is that one cannot obtain the constant 
term in the asymptotic expression, because this would require an exact expression for 
the generating function and an investigation of this exact expression in the vicinity of 
its dominant singularity. Despite this fact, Opheusden en Redig were able to derive 
the constant term in (2.3) by studying the continuum limit of the random walk. They 
conjectured that the asymptotic behavior of the continuum limit is the same as for the 
random walk, and provided strong numerical evidence. Here we provide the proof of this 
conjecture. 

For the voter model in Figure 1(b), our Theorem 1 gives 

f {io , jo) [(X, Y) hits V before %*] ~ i -> oo. 

We next present a result for random walks with a negative drift. 
Theorem 2. Let (X,Y) be a random walk satisfying (H1)-(H4). // 
(2-4) E-i 

there exist constants B(j$) G (0, oo) and p G (0, 1) such that 

■3/2 
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(iojo)i(. X > Y ) llits V before H] ~ B(j )-^, i -)■ oo. 



The same result can be shown to hold for random walks with no transitions to the North, 
North-East and East. Introduce the assumption 

(H2') p_i ) i+P-i,o+P-i,-i+l>o = 1, P-i,i / 0, p± _i ^ andp_i,i+pi _i / 1. 

Note that this assumption is satisfied by the random walk in Figure 1(c), and that (HI) 
and (H2') immediately render (H3) and (H4). We have the following result. 

Theorem 3. Let (X, Y) be a random walk satisfying (HI) and (H2'). There exist constants 
C(jo) G (0, oo) and p G (0, 1) such that 

F M) [(X,Y) hits V before H*] ~ C(j )^, i -»• oo. 
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Figure 2. A tandem queue. 



Explicit expressions for the quantities -B(j'o), C(jo) an d p in Theorem 2 and Theorem 3 
are derived in Section 4. For the model of Godreche et al. [9] in Figure 1(c), Theorem 3 
yields 

(2.5) P(ioj )[(^. Y ) hits V before «1 ~ COo)^, *o -> oo, 
with 

P = j[l-A- V(1-A)(1-2A)]. 
This result matches with Godreche et al. [9, Equation (6.41)]. As in Opheusden and Redig 

[16], Godreche et al. [9] used a functional equation to derive the term p l °i in (2.5), 
but since the functional equation was not solved explicitly, it was impossible to derive the 
constant term C(jo). We provide an exact expression for C(Jq) in Lemma 10. 

Let us finally present a result for random walks with a negative drift in the horizontal 
direction, and zero drift in the vertical direction. 

Theorem 4. Let (X,Y) be a random walk satisfying (H1)-(H4). If 

(2.6) E-i<y<i m,i < and E-i<i,Ki^' = °' 
i/iere exists a constant D G (0, oo) swc/i i/iai 

P (ioJo) [(X, F) hits V before «*] ~ i -> oo. 

An example of such a walk is displayed in Figure 2, with p\ Q = A, po,-i = ^ an d 
p_ 1 1 = i/, assuming A < v. This walk represents the transitions of a tandem queue with 
Poisson arrivals at queue 1 with rate A, exponential services at both queues with mean 
1/u, and all customers traversing from the first queue to the second queue before leaving 
the system. In this case (2.1) describes the probability that, starting with Iq customers in 
queue 1, and jo customers in queue 2, queue 2 empties before queue 1. The constant D is 
identified in Section 4 and for the tandem queue in Figure 2 takes the form 



D 



A 



The remainder of the paper is structured as follows. We first derive, in Section 3, an 
explicit expression for the probability (2.1), for which we rely heavily on earlier work in 
[8, 11]. In fact, this analysis leads to expressions for the generating functions, of which the 
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probabilities in (2.1) are the coefficients, in terms of integrals that involve certain conformal 
gluing functions. In Section 4 we prove Theorems 1-4 by asymptotically evaluating the 
exact integral expressions for the generating functions obtained in Section 3. Finally, 
inspired by [16], we present the continuum limit for the zero-drift random walk in Section 
5. This continuum limit allows for an easy and explicit analysis of the type of event 
described in (2.1). We show that the constant that arises in the continuum limit matches 
with the exact constant we obtain in our precise asymptotic expression for (2.1). This seems 
to suggest an interchange-of-limits, but establishing a formal proof of this fact remains an 
open problem. 

3. Exact integral representations 

We first derive an explicit expression for the probability (2.1) in terms of integral 
representations for the generating functions 

(3.1) h io ' jo (x) = ^> (i0iio) [(X,Y) hits H before V* and at (i,0)]s i-1 , 
/^o jo ( y ) = J^P (ioJo) [(X,y) hits V before U* and at (0, j)]y i_1 . 

m 

Let 

= F {ioJo) [{X,Y) hits (0,0) before U* U V*] 
and note that the probability (2.1) follows from 

(3.2) ¥ {io>jo) [(X,Y) hits V before %*] = h io ' jo (l) + = 1 - /> J °(1). 

The second equality in (3.2) is due to the fact that under (H1)-(H4), the process hits the 
boundary with probability 1, see [7]. In Subsection 3.1 we introduce some classical notions 
regarding the framework in [6] that aims at solving functional equations for the above 
generating functions using the theory of boundary value problems. This framework leads 
to the results presented in Subsection 3.2. 

3.1. Basic properties of the kernel. A common and crucial quantity of interest in the 
study of walks with small steps (as in (HI)) in the quarter plane is the kernel 

( 3 - 3 ) K(x,y) = xy[Y,_i <id ^iPi,jX l yi - 1]. 

It can also be written as 

K(x, y) = a{x)y 2 + b(x)y + c{x) = a(y)x 2 + b(y)x + c(y), 

where 
(3.4) 

a(x) =pi ) ix 2 +p 0) ix+p-i,i, b(x) =pi fi x 2 -x+p-i fi , c(x) =pi _ix 2 +po,-l£+P-i -l, 
a(y) =Pi,W 2 +Pi,oy+Pi,-i, b(y) = p ,iy 2 - y+ po-i, c(y) =p- ltl y 2 +p- lfi y+p- 1 ^ 1 . 
We define 

(3.5) d(x) = b{xf - 4a(x)c(x), d(y) = b(y) 2 - 4a(y)c(y). 

Under (H1)-(H2), the polynomial d has degree three or four, and we denote its roots by 
{^}l<f<4 witn 

1^1 1 ^ \%2\ ^ 1 3^3 1 ^ \x^\, 
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and X4 = 00 if d is a third-degree polynomial. One can easily see that x\ G (— 1, 1) and that 
X4 G (1, 00) U {00} U (—00, —1]. As for the roots X2 and X3, they are positive and such that 
x\ < X2 1 ^ X3. Furthermore, X2 = 1 (resp. X3 = 1) if and only if j^i iPiJ = 

(resp. X)— = 0)- The polynomial d in (3.5) and its roots {yi\i^i^A satisfy similar 
properties. These facts, as well as Lemma 5 below, are proved in [6, Chapter 2]. 

Lemma 5. The polynomial d is positive on (x2,X3) U (2:4, x\) and negative on {x\,X2) U 
(2:3, X4). Similarly, the polynomial d is positive on (2/2,2/3) U (2/4,2/1) and negative on 

(2/1,2/2) u (2/3,2/4)- 

In what follows, we call X(y) and Y(x) the algebraic functions defined by K(X(y),y) = 
and K(x,Y(x)) = 0. With (3.3)-(3.5) we have 



(3.6) X(y)=- m ±r y \ Y(X) --«')±S% 



2a{y) K ' 2a(x) 

where above and throughout, we use the principal determination of the square root, 
defined on C \ (— 00, 0], see [10, Chapter 4]. The functions X(y) and Y(x) both have 
two branches, called Xq, X\ and Yq, Y\. Lemma 5 asserts that they are meromorphic on 
C \ ([2/1,2/2] U [2/3,2/4]) and C \ {[x\, x 2 ] U [x 3 , x 4 ]), respectively 

In the non-zero-drift case, we fix notation letting Xq(1) < Ai(l) and Yq(1) < Y\(l) (see 
[6, Lemma 2.3.4]). Then we have on the whole of C (see [6, Theorem 5.3.3]) 

(3.7) \X (y)\ ^ |Xi(2/)|, |Yo(z)| < |Yi(x)|. 

In the zero-drift case, Xq(1) = -Xi(l) = Yq(1) = Y\(l) = 1, and the previous convention 
has to be changed: we then choose |Xo(— 1)| < |Ai(— 1)| and |Yo( — 1)1 < |Yi( — 1)1 ( see [6, 
Lemma 6.5.1]). Then (3.7) still holds on the whole of C (by a continuity argument). 

Let us finally introduce 

(3.8) M7o(g) = 7 ,v V I /° )d(x) k [-b(x)] j °-( 2k+1 \ 
v ' noK ' [2a(x)]3o \2k + l) v ; L v n 

This quantity appears in the expression of /i l0J0 (x) that we shall give in Theorem 6. It is 
closely related to Y{x): since d{x) is non-positive for x G [xi,^] U [£3,2:4], see Lemma 5, 
the two branches (3.6) of Y{x) are complex conjugates in these intervals. Expression (3.6) 
and some elementary calculations then yield 



(3.9) Y (x) jo - Ytixy = ±2i v / -d(x)/i jo (x). 

In order to determine the sign ± in (3.9), we have to specify whether x —> [x±,X2] U [2:3, £4] 
from above (4,) or below (f): indeed, remember that the branches Yq{x) and Yi(x) are not 
meromorphic on x G [x±, X2] U [2:3, 2:4]. For instance, if x 1 [x\, X2], we have ± = — in (3.9), 
and if x f [2:1,2:2], we have ± = + (see [11, Subsection 3.2]). 

3.2. A crucial conformal mapping. Before deriving an expression for the generating 
function h 10 ' 30 defined in (3.1), we first need to introduce a certain conformal mapping. 
For this, define the closed curve in the complex plane 

X ( [2/1 , 2/2] ) = X ( [2/1 , 2/2] ) U X 1 ( [2/1 , 2/2] ) , 
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which is symmetrical with respect to the real axis (since for y £ [2/1,2/2], Xq(v) and X\(y) 
are complex conjugates, see Lemma 5 and (3.6)) and goes around the segment [xi,X2] (see 
[6, Theorem 5.3.3]). Denote by 

yX([ yi ,y 2 ]) 

the set surrounded by X([yi, 2/2]), which in addition contains [xi,^]. For instance, in the 
case of the simple random walk (with pip = po,i = P-1,0 = P0,-l = 1/4), [2/1,2/2]) 
is the unit circle (see [6, Theorem 5.3.3]), hence <£X([yi, 2/2]) is the unit disc, since 
— 1 ^ xi,X2 ^ 1 (see Subsection 3.1). 

Let us now introduce a conformal gluing function for the set @X([yi, 2/2]), i- e -, a function 
w such that 

(i) io is meromorphic in &X( [2/1,2/2])) an d for t on the boundary of the latter domain, 
\im x -}t w(x) exists (in CU {00}), provided x £ WX( [2/1, 2/2]); 

(ii) establishes a conformal mapping of [2/1, 2/2]) onto the complex plane C cut 
along an interval; 

(iii) For all t on the boundary of &X([yi, 2/2]), i- e -, f° r an * i n -^([2/1,2/2]), wit) = w(t). 

Examples. If the domain @X([yi, 2/2]) has a simple shape (for example, if it is bounded 
by a circle, or by an ellipse), then it is an easy task to determine a suitable conformal 
mapping. 

A first example concerns the case where &X([yi, 2/2]) is the unit disc (this is the case 
for the simple random walk, see above). Then the function 



is a suitable conformal gluing function. Indeed, (i) is obvious; further, w is a conformal 
mapping from the unit disc to C\ (—1/4, 00), so (ii) is satisfied; finally for t = exp(i^) we 
have w(t) = — l/(4sin(</>) 2 ) = w(t), and (iii) follows. 

Another example is when the domain @X( [y%, 2/2]) can be bounded by an ellipse (this is 
the case for the walks, which, as in [9], satisfy (H2')). Then expressions for w can then be 
found in the literature (see, e.g., [6, Theorem 6.3.1]). 

On the existence and the uniqueness of w. The existence (without any explicit expression) 
of functions w satisfying (i), (ii) and (iii) follows from general results on conformal mappings 
(see [13, Chapter 2]). Any linear transformation [aw + j3)/(jw + 5) of w is also a 
suitable conformal mapping. In particular, there is no uniqueness of w. However, there is 
uniqueness of w up to these linear transformations (see [13, Chapter 2]). 

On explicit expressions for w. In most cases, finding an explicit expression for w turns out 
to be challenging. However, for the class of walks at hand, expressions for w are available 
in the literature. In the non-zero-drift case, an expression for w in terms of certain elliptic 
functions is obtained in [11, Section 4]. In the zero-drift case, an expression for w is found 
in [8, Section 2] (it is recalled here in (4.4)), by considering the zero-drift case as the 
limiting case of the non-zero-drift case. We refer to [8, 11, 17] for additional details. 

3.3. Exact hitting probabilities. Denote by id a conformal gluing function as in 
Subsection 3.2. With the notation introduced in Subsection 3.1, we have the following 
result. 
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Theorem 6. For x 6 C \ (x^, X4), 



h l0M (x) = x t0 Y (xy a + - 

7T 



■f2 



w'(t) 



w\t) 



w(t) — w(x) w(t) — w(0) 



y/-d(t)dt. 



The proof of Theorem 6 is given in Subsection 3.4. Theorem 6, together with the fact 
that Yq(1) = 1 (see [6, Equation (5.3.2)]), yields 



(3.10) 



h io ' jo (l) 



1 + 



1 



7T 



■1:2 



w'(t) 



w'(t) 



w(t)-w(l) w(t)-w(0) 



y/-d(t)dt. 



Thanks to (3.2), (3.10) immediately yields an expression for the hitting probability (2.1). 

3.4. Proof of Theorem 6. We only sketch the proof of Theorem 6, because we largely 
mimic the proof in [11] for the case of two positive drifts, i.e., 

A close examination of the proof in [11] makes clear that the result for positive drifts 
remains to hold in the zero-drift and negative-drift cases. The only difference between 
these cases is that the conformal gluing function w introduced in Subsection 3.2 will be 
different. To be somewhat more specific, we now present the four main steps of the proof 
of Theorem 6, closely following the original approaches of [6] and [11]. 



Step 0. Using simple recursion relations it can be shown that h l °' J0 (x), fo i0 ' J0 (y) and 
/iq°q satisfy the functional equation (see [11, Section 2]) 

(3.11) h io ' jo (x) + h io ' jo {y) + - x io y jo = K(x,y)x 



V V>(i 0ijo )[(X(n), Y(ra)) = (X,Y) did not hit U U V between and n] 



Step 1. Thanks to the fundamental identity (3.11), we prove that /i* 0,:,0 (x) satisfies the 
following boundary value problem. Let SfXQyi, 7/2]) be the set introduced in Subsection 
3.2. Then 

(i) h lQ ' 3Q is holomorphic in ^X([yi, 2/2]); 

(ii) For all t on the boundary of SfX([yi, 2/2])) 

h io ' jo (t) - h io ' jo (t) = t io Y (t) jo -t io Y (t) jo . 

The problem of finding a function satisfying (i)-(ii) is a particular instance of a boundary 
value problem with shift (the complex conjugation plays in (ii) the role of the shift), see 
[13] for an extensive treatment of this topic. Items (i) and (ii) follow from [6, Theorem 
6.5.2]. Note that (ii) is easily proved: it essentially suffices to evaluate (3.11) both at Xo(y) 
and X±(y), for all y S [2/1, yi]- In this way, the kernel K(x,y) vanishes, and in fact the 
right-hand side of (3.11) too. Finally, taking the difference of the equations corresponding 
to Xq(u) and X\(y) leads to (ii). 

Step 2. We transform the problem (i)-(ii) into a boundary value problem with a 
boundary condition on a segment. This can be done via a conformal gluing function 
for the set ^X([yi, 2/2]) as discussed in Subsection 3.2. We refer to [17, Section 3] for more 
details. 
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Step 3. The solution of the latter boundary value problem is elementary, see [13] or 
[17, Section 3], and can be formulated in terms of Cauchy integrals. The explicit integral 
representation of h l0 ' J0 (x) follows. A similar expression can be obtained for h l0,30 {y). 

Remark. We can now elaborate on the reasons for assuming (HI) and (H3). First, if we 
would allow larger jumps, we could still obtain a functional equation for the generating 
function of the hitting probabilities, but the technique that is used to solve (3.11) does 
not carry over. Further, finding a conformal gluing function w requires introducing the 
Riemann surface defined by 

{(x,y)€C 2 :K(x,y)=0}. 

For small jumps as in (HI), this Riemann surface has genus or 1, see [6]. If the jumps are 
larger, the genus of this Riemann surface increases, and the problem of finding w becomes 
more intricate. The reason for assuming (H3) is that for pi i i,— l+Pi — l = 1 it is 
possible that the branches (3.6) of the algebraic functions X(y) and Y(x) are meromorphic 
on the whole C. In this case we cannot state (and solve!) a boundary value problem, since 
this requires complex conjugate branches on some interval. 

4. Asymptotic analysis 

We now prove Theorems 1-4, by asymptotically evaluating the integral expressions of 
the generating function h 10,30 derived in Section 3. 

4.1. Proof of Theorem 1 (zero-drift case). Assume (2.2) and let a be defined as in 
(3.4), d as in (3.5), and 

S-i<ij<i y'ftj 



6 = arccos i rv- -2 n/2 -2 11/2 

Lemma 7. The constant A in Theorem 1 is given by 

(4 2) t*M!l 

1 ' ; 2 3 / 2 6>a(l) ' 

In order to prove Theorem 1 and Lemma 7, we first identify an appropriate conformal 
gluing function w for the domain X([yi, J/2])- Define 

d"(x4)/6 + d'(x±)/[t — X4] if X4 7^ 00, 
d"(0)/6 + d"'(0)t/6 if x 4 = 00. 

The next two lemmas are taken from [8, Section 2]. 
Lemma 8. Let f be as in (4.3). In the zero-drift case (2.2), the function w defined by 



(4-3) f{t) 



(4.4) w(t) 




arcsm 



1 2/(t) 



-1/2' 



_3 d"(l) 

is a suitable conformal gluing function for the set @X([yi,y2])- 
This exact expression renders the behavior of w near t = 1. 
Lemma 9. Let w be as in (4.4) and 9 as in (4.1). There exists such that 

a + o(l) 

( 4 - 5 ) w{t) = jY~tW^ *ti- 
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Proof of Theorem 1 and Lemma 7. Thanks to (3.2), it is enough to prove that 

h i <^(l) = l-Aj /i + O(l/i 2 ) 

with A as in (4.2). First, since w(l) = oo (see (4.5)) and since Xi = 1 (see Subsection 3.1), 
(3.10) yields 

(4.6) *<•*(!) _ 1 - i£ e^^-gML.^At. 
For t € [#1,1], we introduce 

(4.7) Ho (t)y^diT) = Y / Mt-l) k , ^JLo) = S A^" 1 )"- 

fc^l w I J fc> _ 1 

The first sum in (4.7) starts at k = 1 because 1 is a double root of d. The second sum in 
(4.7) starts at k = —1 because w has a singularity at 1 of the kind (4.5). With (4.6) we 
then obtain 



1 



(4.8) h t °> 3a (l) = l-- f°[a 1 P-i + (a 2 /3-i + a 1 Po)(t-l) + ...}dt. 

7T J Xl 



Since, for p ^ 0, 



we deduce that 



f l ■ (-l) p v\ 
/ fo ( i_ lfd i = i_^ +0(i 



I 



(4.9) /l^°(l) = 1 - ^i^i + O(*o 2 ). 

TTlQ 

It remains to identify ct\ and First, since w has a singularity of order tt/9 at 1 (see 

(4.5)), it is immediate that /3_i = —tt/9. In addition, since 1 is a double root of d, we 
obtain that «i = fjL jo (l)[-d"(l)/2] 1 / 2 . Moreover, the equality d(l) = together with (3.8) 
implies that 



j [-6(l)]io-i jQ 



£(1) 
o(l) 



Oo-l)/2 



[2a(l)]J'o 2o(l) 

The last identity follows from —6(1) = 2[a(l)c(l)] 1//2 , which indeed holds because d(l) = 
and 6(1) < 0. Moreover, under assumption (2.2) we have a(l) = c(l), in such a way that 
ot\ = [j /(2a(l))] • [-d"(l)/2] 1 / 2 . The proof is completed. □ 

Remark. Theorem 1 provides first-order expansions of /i J0J0 (l) and h %0M {l). By extending 
our approach, we could obtain expansions up to any order, see (4.6)-(4.9). 

Remark. Notice that d"(l) < 0, so that A > 0, see (4.2). Indeed, it is proved in [6] that 
under (H1)-(H4) and (2.2), only two roots of d are equal to 1. In particular, d"(l) ^ 0. 
By continuity of d"(l) with respect to the parameters {pi,j}-i^i,j^i, it is enough to check 
that for one walk, we have d"(l) < 0. This can be easily done, for instance for the simple 
random walk. 
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4.2. Proofs of Theorems 2 4 (negative-drift case). 

Lemma 10. The constants B(j$) and p in Theorem 2 are given by p = X2 and 



3/2 

(4.10) B(j ) - 



2vr 2a{x 2 ) 
with /?o as in (4.14). 



c{x 2 ) 
a(x 2 ) 



(io-i)/2 



d'(x 2 ) 1/2 Po, 



Proof of Theorem, 2 and Lemma 10. We now assume (2.4). For t £ [xi,X2] define 
(4.11) 

With (3.10) we then obtain 



•'.'2 



V°J°(1) = 1 + - / t ,0 [aoA) + (oiA) + aoA)(*-X2) + -.-]v^ r tdt. 

7T 



1 

) = 1 + 

Since, for p ^ 0, 

fX2 3/2+io+p 

(4.12) / f°(x 2 -t) 1 /2 + P dt = r(p + 3/2) 2 3/2+p +O(i - 5/2 - p ), 
we deduce that 

(4.13) h«*(l) = 1 + a0/30r(3/ 3 2 / ) 2 ^ /2+t ° + 0(i~ 

In particular, the fact that p = x 2 6 (0, 1) is a direct consequence of (2.4) and [6, Lemma 
2.3.9]. It remains to identify a® and (3q. First, we obviously have otQ = d'(x 2 ) 1 / 2 pj (x 2 ), 
see (4.11). Further, the equality d(x 2 ) = together with (3.8) gives 



j [-b(x 2 )]^- 1 j 



c(x 2 ) 
a{x 2 ) 



Oo-i)/2 



[2a(x 2 )]i° 2a{x 2 ) 

The analysis of /3o is more elaborate and requires, as can be seen from (4.11), a detailed 
description of the conformal gluing function w for the set &X([yi,y 2 ]). We shall make 
crucial use of a conformal gluing function derived in [6, 11]. Because the description of 
this conformal gluing function would require the introduction of many new symbols, we 
choose to just give some of its most important properties, and we refer to [11, Equation 
(16)] for its full expression. The properties of w we shall use here are the following (see 
[11, Proposition 15]): 

(PI) w has a simple pole at x 2 ; 

(P2) The other possible poles of w are on (£2,2:3) H (X(y 2 ), 00); 
(P3) The set w(Xi([yi,y 2 ])) is a real interval without double points. 
Properties (PI) and (P2) together with (4.11) imply that 

(, 4 - i4 J Po - t: 77-, wTvr 

\\m t ^ X2 [{t - x 2 )w(t)\ 

This quantity is well defined (i.e., finite): thanks to (PI), the denominator of (4.14) is 
non-zero, and thanks to (P2) and since < x 2 and 1 ^ X^{y 2 ) (see Lemma 11), w(0) and 
w(l) are finite. To conclude, let us show that (3q is non-zero. For this, we first note that 
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both and 1 belong to the closure of < &X([yi,y2\). Indeed, X^(yi) ^ 0, see [11, Lemma 
23], and -^(2/2) 1, see Lemma 11. If and 1 are in the open domain < £X([y\,y2\), then 
w(0) 7^ w(l): indeed, by definition (see Subsection 3.2), w is one-to-one in SfX([yi, 2/2])- 11 
either or 1 lies on the boundary of &X([yx,y2]), the latter reasoning still works. If both 
and 1 lie on the boundary and if w(0) = w(l), (P3) gives w(Xi([yi, 2/2])) = RU {00}. 
The latter is a contradiction with the connectedness of @X([yi, 2/2])- D 

Proof of Theorem 4- In this proof we assume (2.6). In particular, this implies that X2 = 1 
and 23 > 1, see [6, Lemma 2.3.9]. Moreover, (PI) then gives that the function w has a 
simple pole at 1, and the identity (4.6) still holds. For t £ [x\, 1], we introduce 

(4.15) ^ o (t) v ^) = ^7^a fc ( i -l) fc , = £ (3 k (t-l] 



k 



Equations (4.6) and (4.15) yield that 



1 



h^(l) = l-- / f°[a /3-i + (ai(3-i + a /3 )(t - 1) + . . .]^==dt. 

TT J X1 s/l - t 

With (4.12), we obtain 

/ i ^(l) = l-^%^i + 0(,- 3 / 2 ). 

Furthermore, the fact that w has a simple pole at 1 implies that = —1, and 

ao = d / (l) 1/ Vio(l)=io^^. 

In particular, the constant D in Theorem 4 is given by D = d' (l) 1 / 2 / (2y/ira(l)). □ 
Lemma 11. The quantity X^(y2) is such that: 

• VYl-i^ij^iiPij < 0> then x e{V2) > 1; 

• ^/E-i^Ki^i = °> i/ien ^(2/2) = 1- 

Proof If X]-i<i j^i = we have 2/2 = 1; see [6, Lemma 2.3.9], and therefore Xg{y2) = 
-2Q(1) = 1, see [6, Lemma 5.3.1]. If we now assume that j^i ipij < and Xg(y2) = 1, 

then thanks to (3.6) we obtain that 2/2 is a solution to 

(4.16) a(y) + b{y) + c(y) = 0. 

On the other hand, it is straightforward that 1 is a solution to (4.16), because this is 
equivalent with the equality Y2—i<i j<iPi,j = 1- Using the root-coefficient relationships, 
we obtain that the other root is c(l)/a(l) 1. This contradicts the fact that 2/2 < 1> which 
is proved in [6, Lemma 2.3.9]. This means that -^(2/2) 7^ 1. To prove that Xg(y2) > 1, 
it is enough to do this for one walk, using the continuity of the different quantities with 
respect to the parameters {Ptj}— i<tj<i- This is easily done for, for example, walks with 
Pl,0 +Po,i +P-1.0 +P0-1 = 1. □ 

Remark. The expression of (3q given in (4.14) can be considerably simplified in some cases: 
see [11, Remark 16] and [17, Theorem 3] for a list of cases where the conformal mapping w 
takes a particular simple expression (those simple cases are particular instances of walks 
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for which a certain group of automorphisms in the sense of Malyshev [15] is finite). As an 
example, for walks with pi ; o +Po,i +P-i,o +Po,-i = lj we have 

(t - xi)(t - x 4 ) 



(4.17) w(t) 
for w 
w(t) 



(t - x 2 )(t - x 3 )' 

As another example, for walks with p—n +Pi,o + Po,-i = 1 as in Figure 2, we have 

t 



(*-X 2 )(t-[p-l 1 lP0,-l/(P? 1 0^)] 1 /2)2- 

Remark. The quantity S(jo) in (4.10) must be positive, since it governs the asymptotic 
behavior of the hitting probabilities, see Theorem 3. As we shall see now, we can actually 
also prove directly the positivity of B(Jq). First, since both (2.4) and (2.6) imply d'(x2) < 0, 
we have that B(Jq) > if and only if /3q < — the latter inequality is not obvious from 
(4.14). In fact, since /3q ^ (see the proof of Theorem 2 and Lemma 10), we can check 
that /3o < for one particular walk, by continuity of (3q with respect to the parameters 
This is easily done with the expression (4.17). 

Proof of Theorem 3. The proof of Theorem 3 is very similar to that of Theorem 2. First, 
X(\yi,y2\) is an ellipse [6, Theorem 6.3.1], and classical expressions for the conformal 
gluing function are available [6, Equation (6.3.1)]. We then transform this mapping by a 
linear transformation, in order to have a pole at x 2 - We call this new function w. Next, 
the integral representation of h 10 ' 30 in Theorem 6 remains to hold (see the four steps of 
Subsection 3.4). In particular, the expansion of h t0,30 {l) given in (4.13) is still true. As 
for the constants, p = x?, /?o i s as m (4-14) and C(jo) has exactly the same expression as 
B(jo) m (4.10). □ 

5. Continuum limit for the zero-drift case 

In this section we derive the continuum limit for the zero-drift random walks that satisfy 
(HI) and (2.2). As it turns out, the limiting process is much easier to analyze than the 
original random walk, and in particular, the hitting probability in (2.1) is obtained in a 
straightforward manner. In what follows, we write f° r j^n % f° r ^+ an d V 

for (hereafter, we shall identify R 2 and C). Moreover, let => denote convergence in 
distribution. 



Theorem 12. 

<t hm+, n ->• oo 



n" 1 / 2 (X(Lntj),y(Lntj) tGR+ (X C ,Y C ) = (X^Y^ 



with the generator of (X c , Y c ) given by 

(r,A) L= \[ f2 "'^ + 2 ^ ijPi ^y + ^ 

Proof. The generator (5.1) simply follows from studying the limit behavior of 

1 



e 2 



[Y,i,jPi,jf( x + ie ,v + j e ) - f( x >y)] 



as e — > 0, for any function / regular enough; see [5], see also [1, Section 5]. Indeed, we 
have the Taylor expansion 

Eij Pi,jf( x + iz,V + je) = f{x, y) + eDf(x, y) + e 2 Lf(x, y) + 0(e 3 ), 
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where L is given by (5.1) and 

D = (5(1) -5(1)) A + (a (i)_ c (i))|.. 

Further, because of our assumption (2.2) on the drift, the operator D is zero. 
Let 
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□ 



(5.2) 



P 



Eij j 2 Pi,j 


1/2 


\a(l) 


+ c(l)" 


1/2 






[5(1) 


+ c(l). 





with a, c, a and c as in (3.4). Investigating the continuum limit leads to the following 
result. 



1 



sm(6)y 
j3x + cos(9)y 



Proposition 13. 

F {x ^ y) [(X c ,Y c ) hits V before H*\ = 
raf/i as in (4.1). 

From Proposition 13 the following asymptotic result can be easily extracted. 
Corollary 14. For fixed y > 0, 

(5.3) F (x v) \(X C , Y c ) hits V before U*) ~ ^ , a; -> oo. 

v p0 a; 

The next result shows that the constant term in (5.3) corresponds with the constant 
term in the asymptotic expression for the original zero-drift random walk. 

Proposition 15. Let A be the constant in Theorem 1. Then 

_ [-d"(l)]V2 _ sin (g) 

2 3 / 2 6>a(l) p0 ' 
Proof. Using (4.1) and (5.2), it is enough to show that 



(5.4) 



[-d"{l)} 1 ' 2 

23/2 



sin(0)a(l) 



5(1) +5(1) 



1/2 



a(l) + c(l). 

where we have used that a(l) = c(l) and 5(1) = 5(1), due to (2.2). Since 

12 \V2 / ,. l2 \l/2 



sin(6>) 



± 1 



taking squares at both sides of (5.4) yields 

-d"(l) = / _ IZijijPiJ? 
8 ^ 4o(l)o(l) 

so that it suffices to prove that 

-d"(l) 



sin(0)[a(l)5(l)] 1 /2, 
l). Since 

Ei,j Wtjf 
4a(l)o(l) 



± 1 



a(l)5(l), 



(5.5) 



MiRi)-EuW.,i]' 



2 —\-/~\-/ u-n,J "Jri,3\ 

Combining (3.4) and (3.5) yields 

d"(l) = 26(1)6"(1) + 26'(1) 2 - 8a'(l)c'(l) - 4a"(l)c(l) - 4a(l)c"(l). 
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Since a(l) = c(l) (see (2.2)) and a(l) + 6(1) + c(l) = (because Y2ijPi,j = 1), we obtain 
d"{\) = 2b'(l) 2 - 8a'(l)c'(l) - 4a(l)[a"(l) + &"(1) + c"(l)]. 

Using that a'(l) + b'(l) +c'(l) = (the fact that the drifts are zero indeed implies that 
the kernel K(x, 1) has a root of order two at 1, see Subsection 3.1), we get 

d"(l) = 2[a'(l) - c'(l)] 2 - 80(1)^1,1 +Pi,o +Pi,-i] = 2[a'(l) - c'(l)] 2 - 8a(l)5(l), 
from which (5.5) is an immediate consequence. □ 
Proof of Proposition 13. As in [1, Section 5] we introduce 
(5.6) 4>(x, y) = a(/3x + cos(0)y, sm(0)y), 

where we have set 

(5 ' 7) a = ¥m^~pp^' 

The motivation of this definition is twofold (for a proof of the facts below, we refer to [1, 
Part 2]): 

• The covariance of (j)(X c , Y c ) is described by the identity matrix, so that (j)(X c , Y c ) 
lies in the domain of attraction of the standard Brownian motion; 

• The function <fi satisfies 

1 d 2 d 2 

(5.8) L(W) = -AW, A = ^2+^2- 

Introduce the notation (U C ,V C ) = (j){X c ,Y C ). Using (5.6) and the usual identification 
between R 2 and C, this random process lives in the cone 

(5.9) {pexp(io;) : ^ p < oo, < u < 9}. 

Another property of the process (U c , V c ) is that (below, exp(i6)T-L stands for the rotation 
by an angle 9 of the horizontal axis ~H) 

(5.10) ¥ {Xty) [(X c , Y c ) hits V before %*] = P^ x , y) [(U c , V c ) hits exp(i9)H before H*], 
and we shall now determine 

(5.11) F (UiV) [(U c ,V c ) hits exp(i9)H before H*]. 

Let us first prove that the probability (5.11) is harmonic inside of the cone (5.9). For this, 
note that for (zq, jo) S 

EijPi^do+iJo+j)^^) hits ^ before W] = F M) [(X,Y) hits V before %*]. 

As a consequence, on the positive quadrant we have 

L(F [X)y) [(X c ,Y c ) hits V before %*]) = 0, 

with L as in (5.1). Thanks to (5.8), the probability (5.11) is then harmonic in the cone 
(5.9). Therefore, we can search for it as the imaginary part of a certain holomorphic 
function H — unique, up to some real additive constants. In other words, 

(5.12) F {u , v) [(U c ,V c ) hits exp(i9)H before %*] = %{H(u + iv)]. 
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To find this function H, we shall exploit the fact that we know its boundary values. 
Since it is impossible (resp. certain) for the process (X C ,Y C ) to be absorbed when starting 
from T~L* (resp. V), one must have 

-»L-«WJ " \ i if zeexp(i0)H*. 

A suitable choice for H is 

(5.13) H(z) = ±]og(z), 

where log refers to the principal determination of the logarithm: this means, see [10], that 
for z £ C \ (— oo, 0] we have 

log(z) = log(|z|) + iarg(z), 

where in the right-hand side, log stands for the classical logarithm function on (0,oo). 
For (5.12) and (5.13) we deduce that 

¥ {UtV) [(U c ,V c ) hits exp(i9)H before %*} = - arctan 

Thanks to (5.6) and (5.10), the proof of Proposition 13 is then completed. □ 
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